Discrete Fourier-based Correlations for Entanglement Detection 
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We introduce two forms of correlations on two d-level (qudit) systems for entanglement detection. 
The correlations can be measured via experimentally tractable two local measurement settings and 
their separable bounds are determined by discrete Fourier-based uncertainty relations. They are 
useful to estimate lower bounds of the Schmidt number in order to clarify generation of a genuine 
qudit entanglement. We also present inseparable conditions for multi-qudit systems associated with 
the qudit stabilizer formalism as another role of the correlations on the inseparability problem. 



Entanglement plays a central role in the development 
of quantum physics and quantum information science. In 
order to understand the entanglement phenomena, it is 
important to ask what is the simplest test for entangle- 
ment. While there are several approaches on this ques- 
tion [lj-Q j an experimentally plausible simple test would 
be formulated from the correlations of the measurement 
outcomes of two local non-commutable observables. A 
reasonable choice of the pair observables is complemen- 
tary ones whose eigenstates form two conjugate bases re- 
lated by a Fourier transform with each other. The corre- 
lation of the canonical position and momentum operators 
is employed for the detection of the continuous- variable 
entanglement 0, H| ■ Similar formalism for discrete sys- 
tems is described in terms of two Pauli operators. 

For the two 2-level (qubit) systems denoted by A and 
B, one of the simplest conditions for the existence of 
entanglement might be 0, H, [f| 0] 



(Z A ®Z B +X A ® X B ) > 1, 



(1) 



where Z := |0) (0|-|1) (1| and X := |0) (1| + |1) (0| are the 
Pauli Z and X operators with the fixed Z-basis {|0), |1)}. 
Associated with the standard quantum key distribution 
protocol, Eq. (JTJ corresponds to the condition that the 
average of the bit error rate and phase error rate is lower 
than 25% Q , namely, it is sufficient to clarify the insepa- 
rability if the average correlation of the Z-basis bits and 
A-basis bits exceeds 75%: 



\ X>'> {i\ + d)Q\®\D (il)> 



(2) 



i=o,i 



where the A-basis is given by the Fourier transform of the 
Z-basis as |0) = (|0) + \l))/y/2 and |1) = (|0) - |1))/V2. 
In order to test this condition, two local measurement 
settings of the measurement in the Z^Ze-basis and the 
measurement in the XaXb -basis are sufficient. Inequal- 
ities (JTJ) and ([2]) suggest that the existence of strong 
correlations both in the complementary aspects can be 
a signature of entanglement. Such types of coexisting 



correlations are also significant for quantification of en- 
tanglement in discrete-variable systems [H as well as in 
continuous- variable systems (Io| . 

There has been growing interest in generating entan- 
glement beyond on the qubit pairs and in verifying that 
the generated states possess genuinely multi-dimensional 
entanglement jllMl4| . To eliminate the possibility of 
lower dimensional entanglement the Schmidt number 



la ] has been utilized in the experiments [12|, |13[ . The 



Schmidt number is defined by the minimum of the max- 
imum Schmidt-rank of the pure states to construct the 
density operator. It indicates the number of degrees of 
freedom contributing to the entanglement and is an en- 
tanglement monotone to quantify the entanglement [l6| . 
An ultimate goal of the experiments is to generate the 
maximally entangled state (MES). The fidelity to the 
MES is called the MES fraction. It represents the close- 
ness of the generated entangled state to the MES. Hence, 
a major theoretical step is to provide accessible tools for 
estimating the fraction of MES and for evaluating the 
entanglement without tomographic effort 



14 



Another interesting topic is detection of entanglement 
for multipartite systems where Eq. ([T]) has already con- 
stituted as an entanglement test in the stabilizer formal- 
ism (3,d,0]- However, the entanglement tests associated 



17, m 



with the multi-dimension stabilizer formalism 
have not been established yet. 

In this Letter we generalize the conditions of Eqs. 
(TTJ) and ([2]) for two d-level (qudit) systems by using 
two different forms of the discrete Fourier-based uncer- 
tainty relations [lH, 0]. We also show that the co- 
existing correlations provide lower bounds of the MES 
fraction and Schmidt number for detection of a genuine 
multi-dimensional entanglement. We further find two 
other inseparable conditions on multipartite systems with 
the stabilizer operators for Greenbergcr-Horne-Zcilingcr 
(GHZ) states and for cluster states on the multi-d-level 
systems. The correlations can be measured by two local 
measurement settings and are thought to be useful in the 
experiments. 

Let us define the generalized Pauli Z and X operators 
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£t V^|j)0'|aiid!:= 



as Z : 

ing a fixed Z-basis {|0) , |1) , • 
conditions \j + d) = \j) and 



1)01, 



7d _ 



E3=, 

■ ■ , |d— 1)} with modulo-e? 
U := 27r/d. They satisfy 
For the qubit 

system (d = 2), the operators Z = and X = X^ 
correspond to the ordinary Pauli operators. For the 
cases of d > 2, Z and X are not self-adjoint but still 
unitary operators. We define the X-basis as a Fourier 
transform of the Z-basis by \k) := Z k \-^= Ej=o li)) = 
Z h \0). The X-basis gives the diagonal expression of 

* = E,td e ~^17)0l- We can verify \(j\k)\ = 1/y/d 
for any j and k. The bases which satisfy this condition 
are called mutually unbiased. 

Discrete Fourier-based uncertainty relations. — In order 
to recall two uncertainty relations [l9|,[2(J associated with 
the Fourier-based relations, let us write the probability 
distributions for Z-basis P(j) := (j\p\j) and for X-basis 
P(j) := (j\p\j). For any two probability distributions we 
have the relation J2j( p2 (j) + p2 U)) < 2 - For the distri- 
butions of two mutually unbiased bases, a strict condition 
holds: dUUU: 



d-l 



E( p2 0) 

3=0 



P 2 U)) < 1 



1 

2' 



(3) 



There is another uncertainty relation concerning the 
complex-valued expectation value of the unitary opera- 
tors 19|, 22 1 . The action of the unitary operator Z (X^) 
applies the phase factor e lu):j onto the j'-th eigenstate |j) 
(|i)), and its expectation value can be seen as the com- 
plex amplitude. The relation (j\Zp\j) = P^e 1 ^ implies 
that (Z) maps the probability distribution of the Z-basis, 
{P(j)}, to a complex number inside the regular d-gon 
inscribed in the unit circle on the complex plane. Simi- 
larly, the relation (X) = P{j) e ~ luJ ^ maps the X-basis 
distribution {P(j)} to a complex number inside the reg- 
ular d-gon inscribed in the unit circle. While the unitary 
condition of Z and X implies \(Z)\ < 1, |(X)| < 1 and 
+ |(X)| < 2, the Fourier-based relation on {P(j)} 
and {P(j)} suggests definite trade-off between the com- 
plex amplitudes (Z) and (X). Intuitively, the amplitudes 
(Z) and (X) cannot be found on the unit circle, simulta- 
neously. A tight condition due to the incompatibility is 
given by Theorem 2 of Ref. [lj| 

| (Z) | cos 6 + | (X) | sin 6 

< i||(Z + Zt)cos0 + (X + Xt)sin0|| =: ||^||, ( 4 ) 

where 9 € [0, n/2] and ||(5|| := max( u | u ) =1 | (u\ 0\u)\. 

Inseparable conditions from the uncertainty 
relations. — We define the MES on a two-qudit 
system as 

d-l d-l 

Vd j=0 Vd j=Q 



and define the total correlation on the Fourier bases as a 
generalization of the operator in Eq. ([2]), 

d-l 

c d E (U) 01 ® li) 01 + 17) (71 ® R) Rl) -(e) 

3=0 

From Eq. (0 it is direct to see that the MES |$o,o) is a 
simultaneous eigenstate of the projective operators: 

d-l 

Edi) 01® li) 0'l)l*o.o) = |*o,o), 



3=0 



d-l 



E(b')01®H>HI)i*°.°> = i $o ^- ( ? ) 

3=0 

This implies that the MES provides perfectly correlated 
Z-basis outcomes and perfectly anti-correlated X-basis 
outcomes, simultaneously. Hence, the central question is 
how high one can enlarge the total correlation without 
entanglement. 

For the product state \(p), we can estimate 

(Cd) 



E,to(nO)^0) + P*ti)PA-j)) 



< 



E^o( p |0) + P*m VE^o(^O') + V(-i)) < 

1 + 1 /d where the first inequality comes from Schwarz 
inequality |a||6| > \a ■ b\. Then, the trivial relation 

Epo'^Vj) = E?=oV(i) and Eq. © lead to 
the final expression. Consequently, for the separable 
state p s = J2iPi \<f>%) 0M ® \ we have (Ai)p 8 = 

E i Pi(C'd)«!'i®v i < E 4 -Pi (1 + 3) = 1 + 2- The equality 
is achieved by, e.g., the product state \(f>) \ip) = |0) |0). 
We thus obtain the inseparable condition for the total 
correlation Cd as c?-dimcnsional counterpart of Eq. @: 

d-l 

Ed-?') 01 ® li) 01 + 17) (7| ® R) R|) > 1 + 5 . (8) 

This condition is equivalent to the condition in Ref. [l4[ 
for two-qudit case with I = 2 and can also be derived 
from the method of Ref. [23j by considering the sum 
ofthe variances Ejl^fli) 0L~ li) 01s) + 52 (|7) <7| A - 
|— i) (— i| B )] with Eq. ([3]) and binding the non-linear 

terms J2j[( p A(j) - Pb{])) 2 + (P A (j) - p bH)) 2 } > 0. 

Similarly we can find another inseparable condition in 
terms of the generalized Pauli operators as follows. From 
Eq. ([5]) it is also direct to see that the MES is a simul- 
taneous eigenstate of the relative amplitudes of the gen- 
eralized Pauli Zs and Xs as (Za — Zb)\&o,o) = and 
(X A - Xjj)|$ ,o) = 0. By multiplying (Z A - Z B )+ and 
(Xa ~ X]j)^ from left side respectively we have the rela- 
tions with the self-adjoint operators, 



2. 



d-l, 



{Z A Zl + Z\Z B )\^ ) 



(X a Xb + X+X^)|$o,o) 



2|*o,o), 
2|*o.o). 



(9) 
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FIG. 1: (Color online), (a) The separable boundary of the 

total amplitude Md = max^ jj^)^ 2 + |{X)^| 2 j as a function 

of the local state-space dimension d. (b) The set of optimal 
distributions that attain the maximum for d = 3, 10, and 20. 



Let us define the total correlation of the amplitudes in 
correspondence with Eq. |T]) as 



Z A Z B 



z\z B 



X A X 



B 



X\X% ) .(10) 



Then, Eq. © implies (Rd)<$ 00 = 2. This is different 
from the situation that \(X)\ = 1 and \(Z)\ = 1 are si- 
multaneously satisfied. Actually, single expectation val- 
ues vanish, (X) = (Z) = 0. It seems impossible to at- 
tain (Rd) = 2 by the separable states although cither 
(Z A Z B + z\Z B ) = 2 or (X A X B + X A X B )_ = 2 can be at- 
tained by the product state, |0)|0) or |0)|0), respectively. 
Hence, the question is how large the expectation value of 
the total amplitude (Rd) can be without entanglement. 

For the product state, we can find a form of the upper 
bound on the total amplitude with Schwarz inequality as 

2(Rd)<f,® v = (Z A )^(Z B ) V + (Z A )^(Z B ) V 

+{x A )^{x B ) v + {x\)^xl) v 

< 2 (\{Z A )^\\{Z B ) V \ + \(X A ) 4> \\(X B ) 



< 2J\(X)^ + \(Z)^J\(X)^ + \(Z^)^ 



< 2 max 

4> 



\(Z) 



\(X)< 



--: 2M d . (11) 



For the separable state, this inequality implies (Rd) p 3 = 

Y,tPi( R d)<t>i® Vt < J2iPi M d = M d- Therefore, we obtain 
another inseparable condition as rf-dimcnsional counter- 
part of Eq. (J): 

- (Z A Z B + Z\Z B + X A X B + X A X B ) > M d . (12) 



We can write M d 

. 2 



max^ 



\(z)^ + \(x)4< 



maxo<9<2L \\Xo\\) by using the uncertainty rela- 
tion of Eq. (01 and the relation ^/|a| 2 + \/3\ 2 — 



m&xg (\a\ cos6 + |/3| sinfl) with a = (Z) and j3 — (X). 
Inequality (fTTj) is saturated by the product of the local 
state that picks up the maximal eigenvalue of the opera- 
tor xe at the optimal value of 9. Figure Q] shows Md from 
d = 2 to d = 100 and three sets of the two Fourier distri- 
butions relevant for the optimal local state. The optimal 
value of 8 is 9 = ir/4 except for d = 3 where Md = 1 is 
achieved at 6 = 0, or 7r/2. Note that we can also obtain 
inseparable conditions with the convex sum of the opera- 
tors p(Z A Z B + Z\Z B ) + (1 -p){X A X B + X\X B ) in that 
case we set a = ^/p(Z) and (3 = \J\ — p(X). It might be 
topical that the two inseparable conditions of Eqs. (|8]) 
and ([T2j) coincide for d = 2 and d = 3 whereas the sets of 
the detectable entangled states include mutually exclu- 
sive subsets for d > 4 as in FIG.[5J An element of the sub- 
sets can be found by considering the noisy MESs V>(p) = 

P$o,o + (1 -P)*Lll,UJ and ^ =-P*o,o + (1 -p)$i,o, 
where we write $z jm = \^i. m ) (®i.m\ with the Bell states 
\$>i, m )AB ■= X A Z B % \$ 0fi ) AB and the floor function |_-J. 
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FIG. 2: (Color online). Noise tolerance threshold p for the 
states 'tpip), 4>(p), and p$o,o + (1 — p)l/d 2 . Circles are for 
the qudit correlations of Eq. ([8| and triangles are for the 
amplitude correlations of Eq. (|12|l . For d > 4, there exist 
parameter regime of p denoted by the symbol X (Y) where 
tp(p) [cj)(p)] is shown to be entangled due to Eq. ^ [Eq. (fl2|) ] 
but its inseparability cannot be verified by Eq. (|12[) [Eq. 



MES fraction and Schmidt number. — We proceed to 
the stronger statement for genuine qudit entanglement 
detection. Interestingly, the two operators take the 
Bell-diagonal forms: C d = Y,t=a *',o + E m ~=o *o,m 
and R d = Efm=o {( cos ^ w + COSTOW ) < ^« 1 }- Using the 
completeness relation Y^i m=o ^i,m — 1; wc can con- 
firm C d = 1 + $o,o ~ E^ m Li < 1 + $0,0 and 

Rd < (1 - COSCj)$ ,o + (1 + COSW) Ef,m=0 *I,m = C 1 ~ 

cosw)$o.o + (1 + coscj)1. We thus obtain two lower 
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bounds of the MES fraction as 
(*o,o) > (Cd)-1, 

(*o,o) > z ((Rd)-(l + COBuj). 

1 — cos u) V / 

Note that the state is uniquely determined to be |$o,o) 
in the case of (Cd) = 2 or in the case of (Rd) = 2. From 
Lemma 1 of [15|, if the fraction is ($o,o) > (k — l)/d, 
the Schmidt number is at least k (See also, Eq. (6) of 
[1H|). Therefore, the Schmidt number is at least k if ei- 
ther (C d ) > l + (fc-l)/dor (Ed) > [(d~k+l)cosuj + (d+ 
k — l)]/d is satisfied. Hence, by setting k — d, we obtain 
the criterion for generation of the genuine qudit entan- 
glement with the Schmidt number of d. One can also cal- 
culate a lower bound of Entanglement of Formation from 
the lower bounds of the MES fraction [25J . Therefore, the 
coexisting complementary correlations have been linked 
to both the detection and quantification of entanglement. 

For the experimental tests, the qudit correlation {Cd) 
of Eq. ((SJ) requires the total number of 2d expectation 
values for the rank-one projections {\j) (j\ <g> |j) (j\} and 
(j| ® |— j) (— On the other hand, the amplitude 
correlation (Rd) of Eq. (fT0|) requires the total number 
of 2d 2 expectation values for the rank-one projections 
{\j) (j\ ® \k) (k\} and Q\ ® \k) (k\}. All elements 
can be obtained by the measurements of the joint Z-basis, 
{\j) \k)}, and the joint X-basis, {\j) \k)}- A standard im- 
plementation could be the combination of the fixed basis 
measurement and the local unitary transformation that 
induces the Fourier transform. It might be instructive 
to emphasize that (Rd) is real-valued by showing the de- 
composition with the projections: Rd = fc{ c °s[w(j — 

*)] \j) 01 ® |fc) (fc| + cosNj + k)\ \j) (j\ ® \k) (k\}.^ 

Entanglement test with stabilizer operators. — Finally, 
we show the conditions for the entanglement detec- 
tion associated with the stabilizer formalism with qu- 
dits [13, 18 1 . Let us define the stabilizer opera- 



T N := xI^Zn [18|, we have (T m _ x + T^_ x ) + (T m + 
Tl) < 2(|(A r t n _ 1 )||(Z m )| + \(Z m _ x )\\(Xl)\) < 2M d for 
the product states as before. This implies another insep- 
arable condition: i|(T m _! + T^_^ + (T m + T^)\ > M d . 
This is a qudit extension of Eq. (12) in Ref. The two 
measurement settings of the TV-joint alternate- Z-and-X 
bases are sufficient for the experimental test. Because of 
the link between conditions of Eq. @ and Ref. 14| , our 
stabilizer method has an advantage to detect the states 
close to 4>(p) over the method of Ref. [HI]. 

In conclusion, we have established a general relation 
between the discrete Fourier-based correlations and en- 
tanglement on two-qudit system. It links two local mea- 
surement settings with the basic concepts on the study of 
entanglement, and the findings could be useful standard 
tools to analyze multi-qudit entanglement. 



tors 



[171, 

for 



181 ] . Let us define the 
X-qudit GHZ state \G) = 



as S\ := 
2,3,-.. ,N 
and St \G) 
ilar to Eq. 



Vd l^k=o \ K ) 
L X fe and S m := Z m _ x Z} n for m = 
The stabilizer conditions S^ \G) = \G) 
\G) hold. For the product states, sim- 



<P>, we have (S x + S{) + (S m + Sp = 

(x N ) + (x\)(xl) ■ ■ ■ + (z m ^)(zl) + 

.„-i)(Z m ) < 2(|(A m _ 1 )||(X m )| + Ki m _i)||(Z m )|) < 
2Md- Then the same bound holds for the full separa- 
ble states, and we have the condition for entanglement 
iK-Si + S\) + (S m + S+,)| > M d - This is the multi- 
partite counterpart of Eq. (fT2]) and corresponds to a 
(i-dimensional extension of Theorem 1 in Ref. (also 
Eq. (11) of Ref. 0). The two measurement settings of 
the TV-joint Z-basis and TV-joint X-basis 0] are sufficient 
to demonstrate the experimental test. Similarly, defining 
the stabilizer operators for TV-qudit cluster state as T\ := 
X\Z 2 , T m := Z m _iX^Z m+ i for m = 2, 3, • • ■ , N-l, and 
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